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Abstract. We show that the black hole solutions of the effective string theory action, where 
%^ I one-loop effects that couple the moduli to gravity via a Gauss-Bonnet term are taken into account, 

f~^ ' admit primary scalar hair. The requirement of absence of naked singularities imposes an upper 

CO . bound on the scalar charges. 
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lyQ. 1. Introduction 

'^ ■ General Relativity describes very well gravity at the classical level, but a quantum 

;_j ' theory of gravity requires the introduction of a more general framework. One of the 

most promising candidates is presently string theory. This theory is believed to change 
drastically the short-range behavior of classical gravity, but also some of its global 
properties can be modified, such as black hole thermodynamics. For example, the study 
of the black hole solutions of effective low-energy theory has shown that, due to the 
presence of non-minimal couplings, non-trivial scalar hair can arise [Q , in contrast with 
classical general relativity, where no-hair theorems [0] rule out this possibility. 

The effects of string theory on gravitational physics are usually investigated by 
means of effective field theory actions, obtained through a perturbative expansion in 
the string tension a. At the tree level, the effective action of the heterotic (but also 
other types of) string contains a coupling of the dilaton with gravity via the Gauss- 
Bonnet term. The black hole solutions of this model have been extensively studied in 
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the literature, both in a perturbative and numerical 0, |^ approach. It turns out that 
the model admits asymptotically fiat black hole solutions with non-trivial dilatonic hair. 
The scalar charge is not an independent parameter, but is a function of the mass of the 
black hole, and is therefore an example of secondary hair [Q . Also the thermodynamics 
is different from that of Schwarzschild black holes. In particular, it was shown that the 
theory predicts a lower bound on the mass of Gauss-Bonnet black holes 0, |[, which 
corresponds to the state of highest (but finite) temperature and lowest entropy. The 
configuration of minimal mass should be identified with the ground state of the Hawking 
evaporation process. 

In order to build realistic models, one should however take into account that in string 
theory other scalar fields are present in the spectrum in addition to the dilaton, as for 
example the moduli, which originate from the compactification of the higher-dimensional 
spacetime. These also couple to gravity through one-loop effects. At leading order the 
coupling term is proportional to the logarithm of a Dedekind r^-function of the moduli, 
which multiplies the Gauss-Bonnet term 0. The effect of the non- minimal coupling of 
the moduli to gravity in a cosmological context has been studied in several papers p[ and 
it has been shown that in some cases it may lead to models without initial singularities, 
but to our knowledge no investigation has been devoted till now to its implications on 
black hole physics. 

On the other hand, it is well known that in effective string actions the 
electromagnetic field exhibits a non-minimal coupling to the dilaton and the moduli 
similar to that of the Gauss-Bonnet term 0. The black hole solutions have been 
thoroughly studied in this case: if one neglects the moduli, one obtains exact 
magnetically charged solutions, with secondary scalar hair, the scalar charge being a 
function of the mass and the magnetic charge [Q. If one instead takes into account 
also one modulus, the general solution can no longer be written in analytic form, but it 
can nevertheless be shown to depend on three parameters [|l^: thus in this case a new 



independent parameter arises, besides the mass and the charge, and one may speak of 
a primary scalar hair. 

In this paper, we investigate if a similar phenomenon can occur in the purely 
gravitational sector. Since we are mainly interested in showing the existence of primary 
scalar hair in the scalar-gravity sector, we consider a simplified model with a dilaton 
and a unique modulus which couples exponentially with the Gauss-Bonnet term. We 
study it both in a perturbative and numerical setting using the techniques developed in 
Ref. and [|], respectively. We find that in fact the qualitative features are similar to 
the case of Maxwell coupling. We also find that an upper limit must be imposed on the 
scalar charges for given mass, in order to avoid naked singularities. This is reminiscent 
of the extremality bounds in the multiscalar Einstein-Maxwell case [^ . 

The structure of our paper is the following. In section 2 we present the perturbative 



solution and discuss its thermodynamical properties. In section 3 we describe the 
numerical solution and the occurrence of an upper bound for the scalar charges. Section 
4 contains a discussion and the main conclusions. 

2. Perturbative Solution 

The bosonic sector of the effective action for the heterotic string in absence of Yang- Mills 
and axionic fields, is given at leading order in oi by 

^— J (fx^^ \n - 2(V$)2 - 2(VS)2 + a(e-2* + 6e-^^)S\ (1) 



where a = a'/8, 5 is a coupling constant of order unity, $ is the dilaton, S is a modulus, 
whose coupling with the Gauss-Bonnet term, S = TZmnpqT^"''^^'' — ^TlmnK'^'^ + Tl^ has 
been taken for simplicity to be of exponential form. The field equations equations can 
be written as 
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and an alougous expression for T^^. The total energy- momentum is conserved but, as 
noticed in [Q for the case of a single scalar, its time component, corresponding to the 
total energy, is not positive definite, due to the contribution of the Gauss-Bonnet term, 
leaving room for the possibility of a violation of the no-hair conjecture. 

We look for spherically symmetric solutions, with scalars $ = $(r), S = S(r). A 
generic spherically symmetric metric can be written as 

ds^ = -A{r)dt^ + ^^dr^ + R\r)dn\ (4) 

Instead of substituting the ansatz (H) into the field equations (g), it is easier to substitute 
it directly into the action, and then vary with respect to the fields <l>, E, A, a and R. 
The action reads: 
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The field equations can then be written in the form: 

{a-^AR^<^y = -2ae-2*[(a-iA'(l - a-^AR'^)]', 

{a-^AR^T.')' = -2a5e-2S[(^-iA'(l - a-^AR'^)]', 

RR" - a-^a'RR' = 2aa-^V'[2a-^AR'R" - (1 - 3a-^AR'^)] 

-2aV"a-\l - a-^AR'^) - R\^'^ + S'^), (6) 

A'RR! + Ai?'2 - (7"2 = -2aV'A'{l - Sa'^AR'^) + AR^{^'^ + S'^), 
A"R + 2A'R' + 2AR" - a" V'(A'i? + 2AR') 

= 4aV'{a-^AA'R')' + AaV"cx~^AA'R' - 2RA{^''^ + S'^), 

where V = (e~^* + (5e~^^) and a prime denotes derivative with respect to r. Only four 
of the previous equations are independent. 

In order to find approximate solutions to the field equations, we adopt the approach 
of Ref. [§ and expand the fields around the background constituted by the Schwarzschild 
metric with vanishing scalar fields, which is of course a solution for a = 0. Our expansion 
is in the parameter a/m?, m being the mass of the background Schwarzschild solution. 
Since a is believed to be of order unity in Planck units, the expansion is valid for large m, 
in the region where aS <C 7^, i. e. for r^ ^ am. For macroscopic black holes (m ^ 1) 
this condition is always satisfied, except in a neighborhood of the singularity, well inside 
the horizon (region of strong curvature). In particular, the approximation is valid for 
the discussion of the asymptotic properties of the fields, and the questions concerning 
the scalar hair. Near the physical singularity, however, the higher order corrections 
to the effective string lagrangian become important and the perturbation theory is no 
longer reliable. We shall however discuss this regime using numerical techniques in the 
next section. 

At this point, it must be observed that the ansatz (^) for the metric is too general 
and still leaves the possibility of a choice of gauge. In order to perform the perturbative 
calculations, the most convenient choice [^ is to impose a = 1, i.e. 

ds^ = -A{r)de + -^dr^ + R\r)dn^. (7) 

This gauge was also used for finding exact charged black hole solutions in effective string 
theory §. 

We expand the fields as follows: 

A = Ao{l + aipi + a'^ip2 + a^i^3 + ■■■), 

R = r + api + a^p2 + a^Ps + • • • , 

$ = a$i + a^$2 + a^$3 + ■ ■ ■ , 

s = j:o + aJ:^ + a^J:2 + a^^3 + . . . , (8) 



where Aq = (1 — 2m/ r). We have normahzed $ such that $ ^ at infinity. This is 
always possible, by rescaling the coupling constant a (this means that our expansion 
is actually in ae"^*"). However, it is not possible to rescale independently also S, and 
hence we take S — > Sq at infinity. The parameters 5 and Sq will always appear in the 
combination Z = 6e~^^". 

Substituting the expansion (||) into the field equations (|^), one obtains at first order 

[r{r - 2m)^[]' = — —, [r{r - 2m)J:[]' = — . (9) 

With the previous boundary conditions , requiring regularity at the horizon r = 2m, 
the scalar fields are uniquely determined at first order: 

Si 1 /l m Am'^\ , , 

$1 = ^ = - + - + TT' 10 

The equations for the metric fields are given at the same order by 

2mi 
p'/ = 0, [(r-2m)V^i]' = -— pi. (11) 

We impose the boundary conditions that p\ -^ const, z/'i — ^ at infinity. We 
are still free to choose the boundary conditions at r = 2mi. Changing the boundary 
conditions at r = 2m yields a reparametrization of the solutions, but no change in 
their physical properties: in particular, the relations between the physical quantities, 
like mass, temperature and entropy, are independent of the parametrization. The most 
convenient choice is to require that the pi and i\}i are regular at r = 2mi. This is equivalent 
to fix the location of the horizon at r = 2m. With these boundary conditions, p\ = 0, 
ipi = 0. This choice greatly simplifies the higher order calculations. 

We can now evaluate the second order corrections. With the stated boundary 
conditions, the equations for the second order perturbations give 

r / X :./!/ 48m^ ^ 
[r{r - 2m)$2] = ^"^i, 

[r(r-2m)S2j = ^ — Si, 

p'^ = -r($f + Sf ) + ^(<|.;' + ZS'/), (12) 

[(r - 2m)^2]' = -(r - 2m)p'^ - 2^-^p'^ - ^p2 
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whose solution is 
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Up to this order, $ and S are proportional in this gauge. However, in order to clarify 
the structure of the solutions, it is useful to go to the next order, even if such corrections 
would of course be modified by taking into account terms of order a^ in the action. A 
long but straightforward calculation gives 

73(1 + ^3) /1\ 73(1 + ^3) /I 
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16480 + 3969 Z^ / 1 \ _ Z (3969 + 16480 Z^) / 1 
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From these results appears that the metric functions are expanded in terms of a^{l + ZY, 
while the scalar fields depend in a more involved way from the parameters. Moreover, 
the functional dependences of the dilaton and the modulus on r are different ii Z ^ 1. 

The perturbative solutions have the following properties: a horizon is present at 
r = 2m, while a singularity is located at the zero of R] the evaluation of this zero 
is however outside the range of validity of our approximation. It can be expected 
nevertheless that for small values of the mass, or great values of Z, the zero can occur 
for r > 2m, leading to the presence of naked singularities. This will be confirmed by 
the numerical results of next section. 

The mass M of the black hole can be deduced from the asymptotic behavior of the 
metric function A and is given by 

M=™ri+i^!(i±^+^^!(i±^v (13) 

V 12 m4 360 m6 y ^ ' 

Its value is greater than that of the Schwarzschild black hole with equal radius. 
Analogously, from the asymptotic behaviour of $ and S one can deduce the scalar 
charges Dq> and Ds, which, in terms of the mass M , turn out to be 

a / 73 a 17110 + 4599^^ a^ \ 
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It is clear that, contrary to the case in which only one scalar field is present 0, the 
scalar charges are no longer function only of the mass of the black hole, but depend also 
on another parameter, which we identify with Z. Hence, in analogy with the dilaton- 



modulus gravity non- minimally coupled to the electromagnetic field ITDI, also in this 



case a primary scalar hair is present in the solution. This gives an example of primary 
scalar hair in pure gravity models. We notice that, at leading order, Z ~ D^^jD^. 

The temperature of the black hole can be defined as usual as the inverse periodicity 
of the time coordinate which renders regular the Euclidean section of the metric. This 
is given by ITT 



T= ^ I ^1 , (15) 

which yields, at order o?: 

r = --J-(l + a2^2(2m)+a3^3(2m)). (16) 

Taking into account (0), a straightforward calculation leads to 

1 / 73 a^(l + Z^) 12511 a^(l + Z3) \ 
SttM V 120 M4 7560 M^ ) ' ^ ' 

The temperature is higher than that of a Schwarzschild black hole of equal mass, but is 
still a monotonically decreasing function of the mass. 

The entropy S can be defined by means of the Euclidean formalism |T^ as 

S = P^-I.. (18) 

where (3 is the inverse temperature and Ie the Euclidean action 

Ie = — / [7^ - 2(V$)2 - 2(VS)2 + a (e'^* + 6e~^'')S)]dV 

16% JM 
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where K is the exterior curvature. A lenghty calculation gives 
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In the range of validity of our approximation, the thermodynamical quantities of 
course do not differ much from their background values, except that now they depend 
on the further parameter Z. They behave differently only for M <^ a, where however 
the approximation breaks down. It is interesting to notice that all the thermodynamical 
quantities are expanded in terms of a^{l + Z^). 

3. Numerical Solution and Naked Singularity 

As discussed previously, the perturbative solution is not valid in the whole domain of 
definition of the solution. Since we are not able to find the exact analytical solution 
of the field equations, we use the numerical approach described in Ref. [ffl. For this 



calculation it is more convenient to choose a gauge in which the metric function R is 
identified with the radial coordinate r, i.e. 

2 

ds^ = -Adt"^ + ^dR^ + R^dn^, 

where A = A(i?), a = a{R). Comparison with (^ yields <j{R) = dr/dR. Of course, 
the physical quantities do not depend on the choice of gauge. However, in order to 
make the comparison of results easier, we give the perturbative expansions in the new 
coordinates: 

A{R) = (1 - 2m/R){l + a^ij2{R) + a'^MR)) + {2m/R^){a^p2{R) + a^PaiR)) + • • 

a{R) = 1 - a^p'^iR) - a^p'^{R) + ■■■ 

$(i?) = $0 + a$i(i?) + aH^iR) + a\^^{R) - ^\{R)p2{R)) + . . . 

E(i?) = So + aSi(i?) + o'T.^iR) + a\T.^{R) - T.[{R)p2{R)) + ... 

where use has been made of the condition ipi = pi = 0, and the functions pi, ipi $j and 
Ej are those obtained above, evaluated at r = i?. In particular, the metric function are, 
up to order a^, 
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In the parametrization ( pi]) the Einstein-Lagrange equations can be written in a 
matrix form (we set all the string coupling constants to be equal to one for simplicity) 

an A" + ai2a' + 0^3$" + cha^" = h, (21) 

where i = 1, . . . , 4 and the entries of the matrices ajj and 6j are 

On = 

ai2 = -a2i? + 4(a'-3A)(e-2*$' + e-2SS') 

ai3 = 4a(A-a2)e-2* 

ai4 = 4a(A-a2)e-2S 

a2i = a^i? + 8Aa(e"2*<l>' + e-2^E') 

022 = -a2(A'i?+2A)-24AA'(e-2*<|.' + e-2^E') 

a23 = 8AAVe-2* 
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62 = -2a%A' + ARi^'^ + S'2)) + 16AAV(e-2*$'2 + g-^^S'^) 
- 8A'V(e-2*$' + e-2Ss') 

63 = 2cr^i?<l>'(A'i? + 2A)-4e-2*A'V 

64 = 2a^i?S'(A'i?+2A)-4e-2SA'V 

The last (constraint) equation is 
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A'{a' - 3A) = 



The numerical integration is performed as follows: we start from a neighborood of 
the horizon and integrate towards infinity. The mass and charges of the solution are then 
evaluated from the asymptotic behaviour of the metric functions and the integration is 
performed again backwards. More technical details on the numerical procedure can be 
found in [||. 

The behaviour of the solution outside the horizon agrees with that obtained by 
perturbative methods in the previous section (see Figures 1-4) and is similar to that 
of the single-scalar solutions discussed in P, |, || . When Dj^ vanishes, we recover the 
solutions described in Ref. [Q. In particular, one must impose a lower bound on the 
black hole mass in order to avoid the occurrence of naked singularities. For non-zero 
value of Ds the situation changes significantly. Taking M and D^ fixed, for small values 
of Dy, the behavior of the solution does not differ much from the case of a single scalar 
field, except that the position of the inner black hole singularity slowly moves up. This 
situation is shown by the solid lines in Figs. 1-4, where the dependence of the functions 
A, a, e~^* and e~^^ against the radial coordinate r is plotted. When Dy. reaches a 
critical value Dj^crit the positions of the singularity Rs and the horizon Rh coincide. 
The dependence of DY,crit against black hole mass is represented on Figure 5. When 
Ds > DY,crit a naked singularity appears. This situation is shown by dashed lines in 
Figs. 1-4. This naked singularity is a continuation of the black hole inner singularity 
and has the same nature. Of course, since the system is symmetric for the exchange 
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of $ and S, an analogous behavior is expected when $ is varied and S held fixed. 
The dependence of DY,crit from the black hole mass M is approximately linear. Hence, 
according to the cosmic censorship conjecture, the mass of the black hole gives an upper 
limit for the modulus/dilatonic field charge. 
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Figure 1. Dependence of the metric function A against the radial coordinate R, Planck unit 
values, when D-^ <^ D-^^^if (thin line) and Dy. 3> D-^^j.^ (thick line). 
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Figure 2. Dependence of the metric function a against the radial coordinate R, Planck unit 
values, when Ds <C Dy,^^^ (thin line) and D^. ^ D^crit (thick line) 



From a mathematical point of view the appearance of this singularity is a 
consequence of the vanishing of the second factor (in brackets) of the determinant Dmain 
of the system (pTI) 



D. 



AA^ + 5A + C 



where 
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Figure 3. Dependence of the dilatonic function e ^* against the radial coordinate i?, Planck 
unit values, when D^ <C Dt^ctU (thin line) and Z^s ^ Dscrit (thick line) 
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Figure 4. Dependence of the modulus function e ^^ against the radial coordinate R when 
I?E ^ DY,crit (thin line) and D^ ^ Dy.^^^ (thick line) 
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A naked singularity occurs when the factor in the bracket of Dmain vanishes before the 
metric function A. Figure 6 represents the dependence of the position of the zero of 
Dmain agaiust the parameter Z. 




Figure 5. The dependence of the critical modulus charge D^crtt against the black hole mass M 
in Planck unit values, for £)$ = 1. 



The thermodynamical parameters can evaluated numerically and compared with the 
perturbative results. This is interesting especially in order tounderstand their behaviour 
for small mass, where the perturbative approach fails. The temperature can be obtained 
from (^) in a straightforward way. For the calculation of the entropy, eq. (|T^ has been 
used. In particular, the Euclidian action Ie was evaluated by adding its definition as 
an additional equation to the main system (PT|). 

The numerical evaluation of the black hole temperature T and entropy S agrees 
with the perturbative results for great M, see the Table 1. A similar agreement holds 
for the entropy. 

Finally, we notice that the thermodynamical parameters stay finite in the extremal 
case. The thermodynamics is therefore analogous to that studied in absence of modulus 
fields 0, 0, except that now one has one further independent parameter (the scalar 
charge) . 
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Figure 6. The dependence of the value of R for which Dmain = 0, against the parameter 
Z = D^/D^. 

Discussion and conclusions 



We have shown perturbatively the occurrence of primary scalar hair in black hole 
solutions of models with more than one scalar field non-minimally coupled to gravity via 
the Gauss-Bonnet term. This result has been checked numerically. From the numerical 
calculations also follows that naked singularities can appear for small values of the mass 
(as in pure dilaton- Gauss-Bonnet models), or for large values of the scalar charges. 
This is a novel feature of the model under study, and can be compared with a similar 
phenomenon occurring in multi-scalar Einstein-Maxwell models [0]. In that case some 
analytical relations for the extremality condition of the black holes can be obtained, 
while in our case this seems not to be possible. We can however conjecture that a 
relation of the same kind exists also in our case. 
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Table 1. The dependence of the black hole temperature T against the mass M and the parameter 
Z « D^/D^ for the numerical and perturbative solutions. The agreement is better for great M 
or small Z, in accordance with the limit of validity of the perturbative approach. 



M 


Z 


T numerical 


T perturbative 


4.0 


0.1 


1.032079 * 10-°2 


9.975080 * 10-03 




1.0 


1.041815 * 10-°2 


1.000250 * 10^02 


5.0 


0.1 


8.094905 * lO-O''^ 


7.966414 * 10-°3 




1.0 


8.169681 * lO-O''^ 


7.974924 * 10-°^ 




10.0 


1.007615 * 10-°2 


9.583720 * 10-°3 


10.0 


0.1 


3.977655 * lO-O''^ 


3.979124 * 10-°3 




1.0 


3.979633 * lO-o^ 


3.979371 * 10-03 




10.0 


4.101194 * 10"°^ 


4.009911 * 10-°3 


20.0 


0.1 


1.986029 * 10-°3 


1.989444 * 10-°3 




1.0 


1.985790 * 10-0-^ 


1.989452 * 10-°3 




10.0 
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